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ON THE HEAT CONTENT OF A POLYGON 


M. VAN DEN BERG * *, K. GITTINS 


Abstract. Let D be a bounded, connected, open set in Euclidean space with 
polygonal boundary. Suppose D has Initial temperature 1 and the complement of D 
has initial temperature 0. We obtain the asymptotic behaviour of the heat content 
of D as time t ^ 0. We then apply this result to compute the heat content of a 
particular fractal polyhedron as t ^ 0. 


1. Introduction 


The conduction of heat or the diffnsion of matter throngh a solid body is of importance 
in the physical and engineering sciences. The classic reference of Carslaw and Jaeger, 
[5], analyses many examples and applications. The mathematical tools used in are 
centred aronnd separation of variables and Laplace transforms and, in many cases, require 
properties of special functions. From a mathematical point of view, the heat equation, 
heat content and heat trace link the underlying geometry of the manifold and its boundary 
and boundary conditions to the spectral resolution of the Laplace operator. Over the 
last few decades, a considerable amount of progress has been made in understanding the 
asymptotic behaviour of the heat content for small time t, see | 12| . 

It was discovered by Preunkert, m, that even in the absence of boundary conditions 
the heat content of a ball B in Euclidean space R"* which is at initial temperature 1, 
while R™ — B has initial temperature 0, has non-trivial asymptotic behaviour as t ^ 0. For 
small t, the initial condition on the complement of B acts in a similar way to a Dirichlet 
0 boundary condition. This was subsequently stated for bounded, open sets with 
boundary in |15| . and proved in |16| . The discussion is simplified by the fact that the heat 
kernel on R"“ is known explicitly. The general situation for the heat content of a compact 
subdomain H in a compact Riemannian manifold M was examined in [3]. The tools 
of pseudo-differential calculus used there rely heavily upon the smoothness assumptions 
on the boundary. Two-sided estimates for the heat content of non-compact sets in R™ 
were obtained in [3]. These estimates are very different from the ones where Dirichlet 0 
boundary conditions are imposed. See [1] and [2]. 

In this paper we denote the fundamental solution of the heat equation on R’" by 


-™/2g-|a!-J/|^/(4t) 


p(x,y;t) = (47rt) 

and for an open set D C R"*, we define 



Then Unix; t) satisfies the heat equation on R"“ 

(2) Aud = X £ R"*, t > 0, 

(see Chapter 2 in m) and 

(3) limu_D(x;t) = 1 d(x), x £ W " — dD , 


where dD is the boundary of D. 
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We define the heat content of D in at t by 

Hoit) = / dxuD{x]t). 

Jd 

So by 0, 

(4) HD{t)= [ dx [ dyp{x,y,t). 

J D J D 

We denote the Lebesgue measure of a measurable set A C R”* by |^|, its perimeter by 
'P{A), and its (m — l)-dimensional Hausdorff measure by (A). 

If D is a bounded, open set in R"*, m > 2, with boundary dD, then Theorem 2.4 
of [TS] implies that 

(5) Ho{t) = \D\ - t ; 0. 

VTT 


In [^, we obtained explicit bounds for Ho (t) for bounded, open sets D in Euclidean space 
with boundary which are uniform in t and in the geometric data of D. These bounds 
imply that 

, 1/2 

(6) Hoit) = \D\ - V{D)— + 0{t), t i 0. 

Vtt 

We observe that if if is a closed set in R"* with \K\ =0 then by the definition of the 
perimeter (see m) and 

(7) \D-K\^\D\, V{D-K)=V{D), Ho-Kit) = Hoit). 

The observations in 0 suggest that 0 holds for all open sets D with hnite perimeter 
and finite Lebesgue measure. The proof of such a statement is well beyond the scope of 
this paper. 

In this paper, we focus on the heat content of a bounded, connected, open set D C R^ 
with polygonal boundary. We introduce some notation and then present the main result: 
Theorem Let 71 ,-, 7 n denote the interior angles of dD. Each such angle 7 j is 
supported by two edges provided 7 j < 27r. By 0 > we may exclude angles 2-^. We label the 
corresponding vertices by Vi, ■ • ■ ,V„ and note that these n vertices need not be pairwise 
disjoint. Let Wj denote the infinite wedge of angle 7 j with vertex Vj such that WjHD 7 ^ 0 
and the boundary of the wedge contains the two edges which are adjacent to Vj and have 
an angle 7 j. Let 

( 8 ) 7 = { 7 i : (sin 7 i)^ < (sin 7 j)^ for all j G { 1 , 2 , • • • , n}}. 

Eor r > 0, we also define the open sector 

(9) Bj{r) = {x&Wj : dix,Vj)<r} 


and 


sup < r : Bi(r) n Bj{r) = 0 for all ^ 7 ^ j, [J Bk{r) C D> . 


( 10 ) 


Theorem 1. Let D C R^ be a bounded, connected, open set with polygonal boundary dD 
with 7 and R as defined in 0 to (|10|| . Then as t fO, 


( 11 ) 


^ 1/2 


Ho{t) = \D\-riD)— + ^5(72)f 


V j = l 


where g : (0, 27r) —>■ R is given 


giP) = 


7 + (1 - I) cot^, e ( 0 , 7 r) U (7r,27r); 
/3 = TT, 


( 12 ) 


0, 
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and k : (0, tt) X (0, 27r) X (0, 27r) 


is given by 


(13) 


k{a, 9,a) = — {—{a + O + a — Tr) cot(a + 9 + a) — {a — tt) cot a) 
2ti 

+ {{a + a — tt) cot((T + a) + (0 + a — tt) cot(0 + a)) , 

In 


for a + 9 + a ^ n, a ^ n, a + a ^ n, 9 + a ^ n, where a denotes the smallest angle 
between We and Wa- In any of the remaining cases, such as a = n, we define k{a,9,a) 
by taking appropriate limits using I’Hopital’s rule. 

The terms which involve the area and perimeter are as expected and agree with those 
in (§. We see that the heat content has a non-trivial dependence on the interior angles 
of the polygonal boundary. 

We observe that /3 i— g{P) is continuous on (0, 2n), decreasing on (0, tt] and symmetric 
with respect to n. That is 

(14) g{l3) = g{2n - /3), 0 < P < 2n. 

By (12 1 and (|14[), we conclude that g is non-negative. We remark that k{a,9,a) is sym¬ 


metric with respect to 9 and a and that k{a, 9, a) = k{2n — 9 — a — a, 9, a). By Lemma 
Section it follows that k is non-negative. 

In addition, if D is a regular n-gon in then 71 = 72 = ■ ■ ■ = 7n = and the 

angular contribution to the heat content is 


n 

-h 2 cot 

TT 


n — 2 


= o( - 

n 


We observe that the coefficient of Laat in the expansion of the heat content of a 
compact domain fl with smooth boundary dH is also equal to 0, see Theorem 1.6 in [4]. 
Here Laa is the trace of the second fundamental form when is oriented with a smooth 
inward-pointing unit normal vector field. 

The expansion for the heat content of a polygon with Dirichlet 0 boundary conditions 
was obtained in [8] . There it was shown that if vd solves the heat equation An = , x € 

D,t > 0 with limtj^o w(x; t) = 1, x £ D and satisfies a Dirichlet boundary condition 
lima,^a;(, n(x; t) = 0 for any xo € dD then 

r +1/2 " 

(15) / dxvD{x-,t) = \D\ — +'S^ c{'yj)t + 0{e~ 

Jd rrf 


-fl2(sin(^/2))2/{32t)s. 


i=l 


where 


c(/3) = 


f 


4sinh((7r — P)x) 


-dx. 


(sinh(7ra;))(cosh(/?a:)) 

We note that both © and ( |15[ ) have angular contributions which are additive. However, 
in Theorem[^there is an additional term in the case where vertices have multiplicity larger 
than 1. No such term is present in (151 since sectors based at the same vertex do not feel 


each other’s presence due to the Dirichlet 0 boundary condition. 

The strategy to prove ( |15[ ) is inspired by im, and relies on some model computations. 
We use an analogous strategy to prove (111. For points x £ D close to a vertex, say 
X £ Bj{r) for some j £ {1,2, ,n}, r > 0, ud is approximated by uw ■ For points 

X £ D which have a distance at least <5 to dD, for some <5 > 0, ud is approximated by 
1. For the remaining points in D, ud is approximated by uh, where H is the half-plane 
which contains D and whose boundary contains the edge of dD nearest to x. As was the 
case in [ 8 ], the model computations involving the infinite wedge Wj are the most difficult 
to carry out. However, in contrast to the case with Dirichlet 0 boundary conditions, we 
must also consider the contribution to the heat content from a vertex which belongs to 
the boundaries of more than one wedge. We deal with these computations in Lemma 
and Lemma [To| in SectionIn Section]^ we carry out the half-plane computations. 

It is quite remarkable that, in contrast to the smooth case , the asymptotic expansion 
in half powers of t of Hoit) in Theorem [^terminates after the term of order t, leaving an 
exponentially small remainder as t 4^ 0- This agrees with the fact that there are no further 
locally computable invariants of D and dD available from which non-trivial quantities 
could be built. A similar phenomenon has been observed for the asymptotic expansion 
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of the heat trace. See for example [7]. The precise form of the exponential remainder 
remains an open problem. The extension of the results in this paper to general polyhedra 
in is another challenge beyond the scope of this paper. 

However, in Section we use Theorem to compute the heat content of a fractal 
polyhedron which is constructed as follows. Let Qo C R^ be an open cube of side- 
length 1. Let 0 < s < 1. Attach a regular open cube Qi^i of side-length s to the centre 
ci,i, i = 1,..., 6, of each face of dQo, and such that all the faces are pairwise-parallel. Now 
proceed by induction. For j — 2,3,, attach N{j) = 6-5^“^ open cubes ,..., Qj,N(j), 
of side-length to the centres of the boundary faces of the cubes Qj-ip,..., (5j_i_jvo-i), 
again with pairwise-parallel faces. We define the fractal polyhedron Da as 


Da = interior 


Qo U 


U U 

i>l l < i < N ( j ) 


for 0 < s < \/2 — 1 (see Figure]^. We note that for this range of s, no cubes in the 
construction of Da overlap. For the two-dimensional construction, see Theorem 4 in [6]. 
In that paper, the critical value s = \f2 — 1, where the squares just touch, is allowed. This 
is due to the fact that the Dirichlet 0 boundary conditions guarantee the independence 
of the heat flow in these touching squares. In this paper, we do not impose Dirichlet 0 
boundary conditions on dDa, and if the cubes touch, then this could give rise to an extra 
term. For this reason, we only allow 0 < s < ^[2 — 1. 



We have that 




1 — 5s^ ’ 


and that the two-dimensional Hausdorff measure of the boundary is given by 


H\dDa) = 6 


1 - 

1 - 5s2 


Both quantities are finite for 0 < s < \/2 — 1. Moreover, the total length of the edges of 
dDa is finite if and only if 0 < s < | and equals 12 In addition to Theorem 1 

in order to compute the heat content of Da, we require the heat content of a sector ot 
angle tt in an infinite wedge of angle ^ where this sector and wedge share one common 
edge and vertex. This will be computed in Section In Section we prove the following 
theorems which give the asymptotic expansion for the heat content of Da as t .1 0. 
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continuous function Ps 


Theorem 2. Let d = | + | • Fix 0 < s < %/2 — 1, s ^ There exists a periodic, 

with period log(s“^) 

„3 


such that 


(16) 




1 + s^ 


- 6 


/ 1 - sM 

1 1 

r i+s\ 

Vl - 5s2^ 

i/F TT 



1 — 5s® 

+ 6pa{logt)t‘^ + (log(t“^))^^^), 11 0. 

It is easy to see that if we write d = , m = 3, in Theorem 


t 


_then d, = 

is the interior Minkowski dimension of the vertices of dDa for 0 < s < whereas for 
I < s < ^[2 — 1, it is the interior Minkowski dimension of the edges of dD^- Below we 
state the corresponding result for the critical case s = |. 

Theorem 3. For s = |, there exists a periodic, continuous function pi 
period log 25 such that, 

(17) 


R with 


132^ 36 

^5 20 5 Stt 57rlog5 


t logt 

+ 6pi (logt)t + (log(t“^))^^^), 11 0. 


In Section below we introduce some further notation and state and prove several 
lemmas. 


2. Additional notation and lemmas. 

Let D be as given in Theorem In the proofs of Lemma to Lemma we use a 
variant of Kac’s principle of not feeling the boundary to reduce the computation of the 
heat content of D to a collection of model computations. For r > 0, <5 > 0, a; G H, A C R^ 
we define 


(18) 

and 


and 


d{x, A) = inf{|a; — z\ ■. z £ A], 

C{5, r) = \x £ D : d(x, dD) < 6, x [J Bk{r) 


D{S, r) = <x £ D : x ^ \J Bk{r), x ^ C{S, r) > . 


Choose R as in (|10[) and let 5 = ^| siny]. 


Lemma 4. 


f dx [ dyp{x,y\t) = \D{5,R)\+0{e ^ 
J D(S,R) JD 


Proof. By [31 Proposition 9(i)], we have that 


1 - 2e 


-s(xy/(st) 


< 


/ dyp{x,y;t) 
J D 


where (5(a;) = min{| 2 ; — y\ : 1 / € R^ — D}. Since x £ D{&, R), 5{x) > 5. Hence 

\D{S, R)\ - 2\D{S, < f dx f dyp{x, y, t) < \D{&, R)\ 

J D(S,R) JD 

as required. 


□ 


We partition the region D — D(^| sin 7 |, i?) into n sectors, Bi{R), of radius R, n 
rectangles, S'^, of height ^|sin 7 | and 2n cusps of height s in 7 |. The contributions to 
Fdoit) from these regions will be computed in Sections |3 4.1 and 4.2 respectively. Each 
sector has two neighbouring cusps. Each cusp is adjacent to a rectangle and a sector. The 
corresponding m-dimensional result to (6) Lemma 7] is the following. 
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Lemma 5. Let D,F,G be non-empty, open subsets ofM.™',m > 2 such that D f] F ^ 
and G <Z D f] F. Let E be a bounded, measurable subset of G. Then 

[ dx [ dyp{x,y,t)= [ dx f dyp{x,y;t)-\-O{e~‘^^^^*^),tl0, 

J E JD J E JF 


where 


Proof. We write 


£ = inf{|a; — y\ ■. x € E,y £ {D U F) Cl dG}. 


f dx [ 

J E J b 


dyp{x,y,t) 


JE JDDF 

(19) = f dx f dy i 

J E J F 

By |T9|, we have that 


and 


y\t)+ [ 

J E 

dx / dyp{x,y,t) 

J DnF^ 

1 

H 

/ dyp{x,y,t) + d 

J E 

J D^r\F 

J E 

> 1 dx 

/ dyp{x,y,t) - 1 

f dx f 

J E J 

E J 

E J b< 

IV 

H 

1 dyp(x,y,t) -2 

F 

'^^^\E\e 

< j dx 

/ dyp{x,y,t) + i 

( dx [ 

J E J 

F J 

E J b( 

< 1 dx 

J E J 

1 dyp{x,y,t) + 2 

F 

”^^^\E\e 


E J DnF'^ 


-p/m 


-p/m 


This completes the proof. 


□ 


Lemma 6. Let i £ {1, ■ • • ,n} and k £N such that i-\-k <n. Suppose 7 ^, 7 ^+!, • • • ,^i+k 
are interior angles of dD which are supported by edges which meet at the same vertex 
Vi = Vi+i = ■ • ■ = Vi+k- Then 
( 20 ) 


[ dx ( 

J BbR) JE 


dx / dyp{x,y,t) 


i+k 


dx f dyp{x,y;t)+ ^ f dx [ di/p(a;, y; t) + 0(e 14 , 

JWj AU., A Jw, 


— I TP. — I \i+k 1 


Proof By Lemma with D = D, F = U’Et^Wy E = U}l^Bj(i?) and G = {z G D : 
d{z, < ^1 sin 7 |}, we have that 

dx / dyp{x,y,t) 


( 21 ) f dx f 

JumBjiR) J D 


We also have that 

( 22 ) 


L 


f dx f dyp{x,y,t) + 0{e ^ t 4 , 0 . 

Ju'mBA(R) J'-l'mwj 

dx / dyp{x,y,t) 

Ju/ffWA 


u]t7Bj{R) 

i + k 


i+k 


= V/ dx dyp{x,y,t)+ V dx dyp{x,y,t) 

j=i-'Bj(R) JWj JWf 

+ 0(e"'^"/'®*>), UO. 


Combining (|2l[) and (22 1 gives (20l. 


□ 
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Lemma 7. Let C denote the half-plane sueh that 0 7 ^ dSjHdD C dH and S-y C H. 
Then 

f dx f dyp{x,y;t)= f dx f dyp{x,y;t)-\-0{e~^ t 4 , 0 . 

j JD J S-y Jh 

Proof. Using Lemma with D = D, F = H,E = S-, and G = {z G D : d{z,Sj) < 
^1 sin 7 |}, the result follows. □ 


Lemma 8 . Let Ci denote the cusp which is adjacent to S~, and Bi{R). Let H be the 
half-plane as in Lemma^ Then 



dyp{x,y,t) 



dyp{x, y-1) + 0 (e «^(=“T')V(32t)^^ ^ ^ 


Proof. Using Lemma with D = D, F = H, E = Ci and G = {z G D : d{z,Ci) < 
^1 sin 7 |}, the result follows. □ 


3. The contribution to the heat content from points close to a vertex. 


In this section we approximate u_d by UWj. We then compute the contribution to the 
heat content of D from a sector, Bj{R) with corresponding angle 77 = /3. We also compute 
the contribution to the heat content of D from two disjoint wedges whose boundaries 
intersect in a vertex of dD. 

Firstly, we define 

(23) Vp{t-,R)=[ dx [ dyp{x,y;t). 

JBjiR) JWj 


Lemma 9. For (3 G (0, tt) U (tt, 27r), 


Vp{t-R) 


fLRff 


2 R ^ 1/2 


+ 





t 


rR\sin0\ . / T \ _\ 2 

— / dx \^—2Rx + arcsin j + R? — x^j e~^ 

+ tiO. 


Proof. Changing to polar coordinates with x = {ri cos , ri sin 61 ) and y = (r 2 cos 62 , T 2 sin O 2 ) 
in p3| ) gives that 

Vp{t-, R) 

r0 /’/3 rR poo . ^ 

= (47rt)“^ / dOi / d 62 / dri / 

Jo Jo Jo Jo 

where A = cos(0i — 62 ). The change of variable r 2 — riA = p gives that 

Vp{t-, R) = {4:Tvt)~^ f d9i f d02 [ dry f dr 2 {r-ir 2 )e~^’'^~^^"^^ /(4t)->'i(i-A )/(4t) 

Jo Jo Jo Jo 

= (47rt)~^ f d9i f d02 f r dr f dp{p + Ar)e~^ /(4i)-’’ (i-a )/(4t) 

Jo Jo Jo J-Ar 

= h+l2. 

We have that 

h = (47rt)"^ d9i d 92 [ rdr [ dppe~'’^ 

Jo Jo Jo J-Ar 

( 24 ) = 

TT 
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and 


rP rP fR /-oo 

I 2 = (47rt)~^ / d9i / d02 / dr Ar^ / dpe 
Jo Jo Jo J-Ar 


-p2/(4t)_r2(l-A2)/(4t) 


P rl 

d9i 

' d02 [ 

R 00 

drAr^e~^ )/(4t) / 

dpe ” 

Jo 

Jo 

Jo 


f d9t 

r de2 


f dp 

Jo - 

’0 

Jo J 

0 

P rl 

dOi 

de2 [ 

R 00 

drAr^e~^ )/(4t) / 

dpe-” 

Jo 

Jo 

Jo 


f d9^ 

r de2 



JO ^ 

'0 

Jo J 

0 

f d9^ 

r de2 



JO ^ 

•0 

Jo J 

r 


= I 3 + li + h- 

Via the change of variables 9i — O 2 = —rj and integrating by parts with respect to 9, we 
obtain 


2g-rVl-A^)/(4t) 


rP nP nR 

/ d9\ / d02 / dr Ar 

Jo Jo Jo 

fR 

(47ri)“^^^ / / c /02 / rfr cos(0i — ^ 2 ) 

7o 7o 

fP—^l fR 2/ 

(47rt)~^^^ / / drj dr cos r/e~^ 

7o ^0 ^0 

(47ri)~^^^ f dO f dr] f dr cost] e 
Jo Jo Jo 


f 


dpe-” 

-r^(sin{ei-02)^)/(4d 


— (sin 77) /( 4 t) 
(sin 77)^/( 4 t) 


(25) 

= ,fl(47rt) J 

dr] / dr r^ cos r]e~^ 
0 Jo 

f0 fR 



/ dOO dr r^ cos 6 e~ 

'0 Jo 

Similarly, 

(26) 

f0 fR 

Ia = P(A7vt)~^^^ / drj drr^\t 

Jo Jo 


1 jn 

— {4:7vt)~ ’ / d6 9 drr' 

Jo Jo 


(sin77) /(4t) 

51 ^-r2(sine)V{4d 


We first compute I 3 + I 4 , and then we deal with /s. By (251 and (261, we have that 

rR 


(27) 


73 + 74 = 2(47rt) f drr^ f dr] {jj — rf) cos t] e 
Jo J B 


— r*^ (sin Tj)'^ / {At) 


where B = [0,/3] n ([0, U [^,27r]). First suppose fi € [0, f], then by (27l we have that 


/3+/4 

= 2(47rt) 


rR r 

' / drr^ j 
Jo Jo 

rR r\ 

= 2(47ri)~^^^ / drr^ 

Jo Jo 

‘ [ drr^ [ 
Jo Jo 


— r^ (sin rj)^ /{At) 

dtp {P — arcsin^) e”"" 


dp {P — p) cos p e 

Sin 0 


= 2/3(47rt) 


' „ „ rR rsin0 

d'il)e~^ ^ — 2(47ri)~^^^ / drr^ dijjilje 

Jo Jo 


-r^-iir / {At) 


rR rsin p 

2(47rt)~^^^ / drr^ / dtp {arcsintp — tp)( 

Jo Jo 
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fH poo pH poo 

= 2/3(47ri)~^^^ / drr^ / d'ipe~^ — 2/3(47rt)~^^^ / drr^ / di/)e 

io io Jo isin/3 




+ 4 i( 47 r/:) / dr (e 

Jo 

1/9 




-r'’(sin/3)'’/(4t) 


- 1 ) 
H 


r2v.V{4t) 


^sin p 

— 2(47rt)~^^^ / dip (aicsinIp — Ip) / drr^e~ 

= ^-2/3(4..)-/^ r d^ 

2 y.in /3 Jo 04 sin/3 

ps\n ^ poo 

— 2(47rt)~^^^ / dip [aicsinip — Ip) / drr^e~'^ 

Jo Jo 

+ 2(47rt)~^^^ f dip {aicsinip — Ip) f drr^e~^ (=>1/3) /(s*)^ 

Jo Jfi 

_ ^ _ ^1/2 
2 0F 


= ^-^e'^-2ptf :^ + ^j-2ti dip 

J si 


(pp ^ 2 ^ r‘"^ arcsini/; - ip 

p p^ sin ^ Jo 


V V'® 


^sin /3 


+ 2(47rt) / dp (aicsinp — p) f drr^e ’’ /(s*)^ 

/o Jfi 


(28) 

FtR^ 2 R I ^sin /3 ^oo 

=-+ (cot/3)i: + 2 ( 47 ri)~^'^ / di//(arcsin?/) —')/)) / drr' 

2 V^r Jo Jfl 

+ t;o. 


2g-r=^/AV(4t) 


Similarly, for /3 £ 


we have that 


/3+/4 


- 1/2 


I drr dri {pi — rp cos rj e 
0 Jo 


— r^(sin77) /(4t) 
r2v>2/(4t) 


-r'*i/’^/(4t) 


= 2(47rt) 

/“i? /“I 

= 2 ( 47 ri)~^^^ / / c/'?/'(/3 — arcsin' 0 ) e 

io ^0 

z' /'I Z' f ^ 

= 2y3(47ri)~^^^ / drr^ / dipe~^ ^ — 2(47rt)~^^^ / drr^ / 

Jo Jo Jo Jo 

p H p 1 

— 2(47rt)~^^^ / drr^ / dp { arca'mp — p ) e ~^ 

Jo Jo 

= 2/3(47rt)~^^^ f drr ^ f dpe ~'^ — 2/3(47rt)~^^^ f drr ^ f dp e ~^ 

Jo Jo Jo Ji 

-I-4t(47rt)~^^^ f dr { e ~^ — 1) 

Jo 

— 2 { 4 ' Kt )~^^^ [ dp {accsmp — p ) f drr ^ e ~^ ^ 

Jo Jo 


/3i?2 


dp 2R 1/2 


-2Pt j 22L _ + 2t-2t dp 

Jo 


2 ""'"./l V’® 04 


arcsin p — p 

p^ 


-I- 2 ( 47 rt) f dp {axcsinp — p) f drr^e ^^"^*^-1-0(46 /(**)) 

Jo Jr 


(29) 

_ pR^ ^ 2R ^1/2 j_ ( p 

2 

02/-^- o',2 


/ TT \ /“I Z*'^ 

^ 1/2 _|_ ^_ _ J _l_ 2(47r/;)“4/2 / dp {arcsinp — p) / drr 
-f UO. 


2„-r-2v>2/(4t) 
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In order to have (291 in the same format as (281, we replaced 1 by |sin/?| in (291. For 
P G [^, 27r], we have that 

I3 + I4, = f drr^ f drj {P — rj) cos rj e~'" 

Jo Jo 

rP 


rH rp 

+ 2{4'Kt)~^^^ J drr^ j drj {P — r]) cos rj t 


-r^(sinr7)^/{4t) 


— Jl + Jo- 

Now by ( |29[ ), 

(30) 

Jl = - + (— — p') t + 2(47rt)~^^^ f dp) (arcsini/) — tp) f drr^e~^ 

2 ^TT \2 J Jo 

+ tiO. 

Via the change of variables r] =2 tx — rj, we also have that 

j,2 (Bin rip/(4t) 


rP 

J 2 = 2(47rt)“^^^ y drj {P — ij) cos rj J dr 

= 2(47rt)~^^^ f drj {P + rj — 2 tt) cos rj f dr 

Jott-p Jo 

t f drj (P + rj - 2Tr) 

Jorr-p (smr?)3 


^2 ^-r-‘(ainri)‘‘/(4t) 


= 2t 


(31) 


— 2(47rt) / dr] {Pr) — 27 v) COST] I 

727r-/3 Jn 

= - ^-cot,5^ t + tiO. 

Hence by (30l and (31|, we see that for P G [^,27r] and t 0, 
(32) 

h+l 4 = ^- + (27r - 2/3 - cot P)t 

2 y/TV 

r-l sin/3| 


dr g-^^CsinJ7)2/(4£) 


r\sinp\ POO 00 

+ 2(47rt)~^^^ / (i?/’(arcsin'i/j —'0) / drr^e~^ ^ /{8t)^_ 

^0 Jr 


to compute I5. Via the change of variables p = rp', we see that 

-(4711) ^ 


1 

d 6 »i 

f 

t 

d02 

r dr r dp 

J 

0 





-(4711)“^ 

(‘ 

1 

d 6 »i 

f 

1 

d02 

r dr a 2 ^ 3 ^-. 2 (i_a 2 )/( 47 ) r ^-PA%-/(4t) 

J 

0 




Jo Jl 


/■I 

1 


1 


-(8711)“^ 

f 

d 6 ii 

/ 

dOo 

/ drA\ 

J 

0 


/o 


^0 ^1 

-(8711)“^ 


ddi 

r 

dOo 

f dp f 

J 

0 


/o 


7o 

+ ( 871 !)“^ 

r 

dOi 

f 

dOo 

f dp [ 

J 

'0 


Jo 


Jl Jr 2 

f'dS 

TT Jo 

'1 

f de2 
'0 

( “"(.IV AO-+0(«-"GO. 


(33) 

By 2.173 (1) and 2.172 in [13], we see that 
dp —1 


1° 


+ 


arctan 


/Vl -242 


(A2p2 + 1_A2)2 2(1 -H2) 2|^|(1- A2)3/2 |A| 
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Therefore by (331, (341 and the change of variables O 2 — 0i = a, we obtain that 

^2 


rP rP coo 

Jo "'Jo */. '‘f‘-AV + l-A‘Y 

It , f , n2 cos a . 

— I do I da \ — (cot(T) + ^ arctanitan(j) . 

Jo Jo V (smcr)3 J 


Thus 


h = - 


— / d6 f da ( —(cot a)^ + arctan(tan cr)^ + 0(te ^ i 4^ 0. 

^ Jo Jo V (smo-)'^ J 


We note that 


where 


arctan(tan(7) = a + U{a), 
'O, ifaG(0,f); 


Hence it is necessary to compute 

-^-ifdefdaUcotaf 


If;3G (f,f),then 



if G 


[-271, 

if G 

G (f ,271). 

G COS G 

\ . 

-?r-( 

(sinCT)s ^ 

) 


6 


2(sin6l)^ 


2t 

TT 


= ( _ ^cotp+- ) t. 

TT TT TT , 


de [ daU{a)j^^^ =2t f dO f dcr-^^^ = (^ + cot ^ t 
Jo Jo (sina)3 (sina)3 V 2/ 


and if /? G (^,27r), then 


2t 

TT 

= 2t 


Cl 

/ de / daU(a) — 
Jo Jo (si 


in (t)s 


[ ^ de f da +2t f dO [ 

J^ Js (sincr)3 J J 


(sin (t)s 

= (2/3 + 2 cot /3 — 2'K)t. 


. cos a 
da — -TTT + 4t 


(sin (j)s 


I>L 


da- 


3ir J 371 (sincr) 


Hence 


(35) Is = 


- iicot/3+ i t 

TT TT r- — » 


/3+i)^, if/3G(0,f); 

-^ + (l-^)cot/3 + /3-f+ l)t, if/3e(f,f); 

+ (l - ^) cot^ + cot/3 + 2/3 - 271 + t, if /3 G (^,27r); 


Therefore, by combining ( |24[ ), ( |35[ ) and ( |28[ ), ( |29[ ), ( |32[ ) respectively, we obtain that, as 

tiO, 


+ 2(47rf) 


- 1/2 f 

Jo 


I sin/3| 


d'lp (arcsin ip — ■ip) / drre 


r 


+ 0{te 


— R? (sin ^)^/(8t) 
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Via the change of variables p = rip and integrating by parts with respect to ip, we have 
rl ain/3| 


2 ( 47 rt)~^^^ / dip {arcsinip — Ip) I drr^e~’^ 
Jo Jr 


= 2(47rt)' 


- 1/2 c 

Jo 


arcsini/) — 1 /)\ f 

J Jr 


dip 

"I sm/3| 


V 

r I am PI 

= —(47rt)~^'^ii® / dip {arcsinip — ip)e 
Jo 

dip ( 1 - - l/>2 


dp p^ e 


Rp 

-R^p^Kri) 


(47rt) f 

Jo 


(36) 


\ y'l - ■!/>2 

'•I sin^l 


dpp^e ^ + 0{te 


-R^{sinpy/{8t)\ 


' Rp 


r \ sin Pi 

= —(47rt)~^^^ii® / dip {arcsinip — Ip)e 

Jo 

A sinPI / 1 

+ (47rt) ^ / di/) I - ^- 

Jo \ pi - ip'^ 1 + -y/l - 


-R^p'^Urt) 

1 


POO 

/ ' 

J Rp 


dp p e 


2„-pV{41) 


Again integrating by parts with respect to ip, we see that 


(47rt) 


- 1/2 


I 


I sin 0\ 


(37) 


dip 


/"I sin , 0 1 _ ^ „ 00 

= (47rt)“^/^ii® / dip {-iPP1-iP^ + iP)e~^ t i 0. 

Jo 


pi — Ip "^ 1 + pi — Ip "^ 


POO 

/ ' 

J Rp 


dpp^e ^ 


Hence by (361, (371 and the change of variables ip = J), we have that 

'f 

Jr 


2{4'Kt) J dip {arcsinip — Ip) j drr^e '' 

^H| sin 


I sin 0\ 


dx (i?^arcsin(|)-fe)e-^^/(^*) 


/•K|sin/3| . _ \ o o o 

+ ( 47 rt )"^/2 J dx (^-xPR'2 - ®2 + Rx^ e~^ + 0{te~^ 


^i?|sin/3| . / rp \ _\ 

(38) = —(dirt)”^'^^ j dx y—2Rx + arcsin j + xpR? — x'^j e 

+ no. 

This completes the proof of Lemma 


I (it) 


□ 


In Section 4.3 we deal with the remaining integral in (381. 


Lemma 10. Let VKi,W 2 be two disjoint wedges in with eorresponding angles 71,72 
respectively such that dWi n dW 2 = {Vi} for some i £ {1,..., n}. Let a denote the angle 
between Wi and W 2 such that 0 < a < tt . Then 




dyp{x,y,t) = fc(a, 7 i, 72 )t, 


where fc(a, 7 i, 72 ) is as defined in (13l. 
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Proof. By changing coordinates to polar coordinates, as in the proof of Lernmaj^ we have 
that 


[ dx [ 

Jwi Jw2 


dyp{x,y,t) 


Af 

= - dOi 

rr ./n 


'71+0 
^71+q: + 72 


dOi 


/‘71 /’71+Q:+72 o o 

(dTTt)-^ / d0i / dOi / dri / 

Jq J'y^+ot Jo Jo 

dr2{nr2) cosie^-e^) 

^ppg-of(l+P^-2p 003(01-02)) 

J, g-oil+P^-Sp 008(01-92)) 

P 


rii r 

/ d0i 

Jo J-1 

rii r 

/ d0i / 

./ 0 «/ 'v 


71+q: 
71+q: + 72 


d92 


71+q: 

71+q:+72 


d6ic 


^OO /’O 

/ / 

70 ^0 

/*oo /’O 

/ f/ri r? 

Jo Jo 

POO PC 

/ rfpp / 

7o ^0 


71+q: 

/’71 ^7l+'^ + 72 

/ dPi / dP2 / dp' 

Jo J -\-ot 


'71+0 

r71 C71+o^ + 72—f^l C 0 )O 

I ddi I dp I dp 

Jo J ^i+a — Oi Jo 


0 (1 + p 2 - 2 pcos( 6 ii - 02 ))^ 


/ 71 +Q —01 
f <7 + 72 


0 (1 + p2 — 2pcosp)2 


7172 

TV 


^71+Q r<T+72 ^00 

I da I dp I dp' 

J cx. J a Jo 

p^l+OL ^o-+72 poo 

I da I dp I dp 

J a J a Jo 

91 r71+0! /■ct+72 o 

f H-/ da dp j 

TT Jo Jo J- 


(1 + p^ — 2 pcos 


+ 


(1 + p 2 — 2 pcos (1 + p 2 — 2 pcosp )2 


dp- 


cos p 


7172 ^ ^ ^ 

TT TT 


r7i+Q ro-+72 

«/ Q J a 

^71+Q ^ 

«/ Q «/ CT 

^71+Q ^ 

«/ Q «/ <T 

y-71+Q ^ 

y 

«/ Q «/ cr 


+ 


2 t 


+ 


2 t 


-C 08 V> (p 2 + (sinp)2)2 

dp 


, cos QP 
d(pT^ 


f 




sinpl^ Jq (p2 + 1)2 


TV TV 


sin Lp\^ 

< 7+72 


/ 


(p2 + 1)2 


, cos ip 
dp-^ 


r 


^(7+72 


<+ 77 ) 

I sin(/p|'^ 


sin 

arctan(| cot (^|) 


dO (cos 6) 


f 

Jo 


do (cos 6)^ 


(39) 


7172 ^ ^ ^ 

TV TV 


^7l+Q ^<7+72 

y “'^y 

«/ Q «/ <T 


^ TT cos (/p 
4 |sinp|3 


+ — 
TV 


2t ni+a PO -+72 /sin (2 arctan(| cot pD) arctan(| cot p|) \ |cosp| 


071+0 o 

y 

«/ Q «/ cr 


dp 


We see that 


Hence 


S{p) := arctan(| cot p|) = 


+ 


|p-f|, if 0 < p < tt; 


I sinp|3 ■ 


|p- 


if TT < p < 27r. 


sin( 2 arctan(| cot pD) |cosp| 1 /cosp 
4 I sin p|3 2 \ sin p 


Therefore (391 becomes 


(40) 


7172 ^ ^ 2t 

TV TV 


^7l+Q PO 

y 

J cx. 'J a 


'^"'"'’'2 /TT cosp 1 f cosp\^ 5'(p)|cosp| 




4 I sin p|3 2 \ sin p 


2 | sin p|3 


By considering the function 


TV cosp 1 /cosp\^ S'(p)|cosp| 


+ d 


4 I sin p|3 2 \ sin p 


2 | sin p|3 
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for Lp in each of the four quadrants, we have that 


2 /, y (sinv5)3 (sinyj j j 


is a primitive for f{^) for all (p € (0, 27r). Integrating by parts with respect to we obtain 


^ r<7+y2 / 

-J, 


/ X cos p 

(tt - (p)-r-. -^ + 

(sin py 


cos p 
sin p 


(41) 


1 / (72 + (T — 7 r) , . (cr — 7 r) „ i 

7 -;-177 “ COt (72 + cr) - —-^ + cot ct - 272 . 

4 V(sm (72 + cr))2 (sincr)2 ' 


Hence, by (401, (411 and integrating by parts with respect to cr, we have that 


dx / dyp{x,y\t) 
JWo 


JW-i JW 2 

= ^ (—(72 + 7 i + a — tt) cot (72 + 7 i + a) — (a — tt) cot a) 

ZTT 

+ 7 ^ ((72 + a - tt) cot (72 + a) + (71 + a - tt) cot( 7 i + a)) 

ZTT 


as required. 


□ 


4. The contribution to the heat content from points close to an edge. 

In this section, we consider points x G C{^ \ siny], R). We partition this region into n 
rectangles S-^ and 2n cusps Ci, each of height sin 7 |. We approximate ud by uh, where 
Tf C is the half-plane such that 0 7 ^ dS-y n dD C dH and S~, C H, as in Lemma 

4.1. The contribution to the heat content from a rectangle. We first compute the 
contribution to the heat content of D from a rectangle, Sj, of height siny] and length 
L, where L G R, L > 0. We have that 

f dyp{x,y;t) = f dy 2 f dyi{4TTt)~^e~^^^~^^^ /(4t)-(a:2-H2) /(4t) 

Jh 4-00 Jo 

= 1 - 

J Xl 


Let Xl G (0, ^1 sinyl) and 2:2 G (0,L). Then, by integrating by parts with respect to xi, 
we obtain 


dx dyp{x,y-t) 
J Sy J H 

[ dx r 

J S-y J Xl 


— I S'-, I 


dC ( 47 rt) ' e 


-l/2^-CV{4t) 


rL, n ^ I sin 7I noo 

|5't,| — (47rt)“^^^ / dx 2 / dxi / 

Jo Jo J Xl 

Jo Jo J Xl 


d( e 


-CV(4t) 


(42) = - ^ 7 ^ 1/2 Q^^l/2^-RREinjf/(32t)-^^ ^ ^ 


4.2. The contribution to the heat content from a cusp. We now compute the 
contribution to the heat content of D from a cusp. Let Ci denote the cusp which is 
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adjacent to S-y and Then 

’§|sin7| 


r r p^\sin-f\ - roo 

/ dx dyp(x^y;t) = / dx {R — y — x"^) / dy (Ant) 

Jci Jh Jo Jo 

^f|sin7| ^- / POO 

— j dx{R— ^R? — x‘^) 1^1 — y dC (47rt) 


-l/2^-|a;-yr/(4t) 


-l/2g-CV(4t) 


p ^ I sin 7 I POO 

= \C,\- dx{R- -x'^) / dC 

Jo Jx 


/m 


Integrating by parts with respect to x, we obtain 

r-^|sin7| _ /lOO 

dx{R — B? — x^) / d^" (47rt)~^^^e~^ 

J X 

"#|sin7l / d2 


2,7(44) 


- 2 , 7 ( 44 ) 


, — (47rt)“^^^ J dx ^Rx — ^ arcsin — ^\/_R2 _ ^2^ e“ 

_^0(g-fl7=in7)7(324))^ ^^0. 


Hence the contribution from each cusp is 

dx dyp{x,y,t) 

JCi J H 

, rf^l sin7l / p 2 ^_\ 

= |7i| + ( 47 rt)“ / dx i —Rx + — arcsin ^ —j + -\/R'^ — x^j e 

_^ 0 (g-fl={sin 7)7(324))^ U 0. 

4.3. The terms from the cusp and sector contributions. Finally, we deal 

with the remaining integrals from the sector and cusp contributions. Each sector has two 
neighbouring cusps so we are interested in the following integral 

2 ( 47 rt)~ ^ j dx ^—Rx H —— arcsin yyyj + — — x^j e 

p I sin 71 

(43) = ( 47 rt)~^^^ J dx ^—2Rx + R^ arcsin + 4 C\/R^ — x^^ e 

By definition of 7 , in ([^, we can write the remaining integral from the sector contribution, 
p^ , as 

T/?|sin/3| . /T\ _\ 

— ( 47 rt)~^^^ J dx y—2Rx + R^ arcsin yyyj + a;\/ R^ — x^j e~ 

( 44 ) = —(47rt)~^^^ / dx y—2Rx + R? arcsin j + R^ — x^j e~ 

_^Q(g-fl3(.in7)2/(8t)), ^^ 0 . 

Adding and ( [44| ), we obtain 

p ^ I sin 7 I ^ / 'V \ _\ 

(47rt)“^^^ J dx \^—2Rx + R? arcsin yyp + x\/R? — x^J e~ 

pR\ sin 7 | , / ry. \ _\ 

— (47rt)“^^^ J dx y—2Rx + R? arcsin ^ —j + x\/R? — x^j e~ 

/ •■H|sin7| , / rp \ _\ 

dx (—2Rx + R^ arcsin ( — ) + x\/R?^-^) e~ 

f|sm7| ^ ' 

^(7(e-«7Bin7)7(324)), 


- 2 , 7 ( 44 ) 


7(44) 


0-2:7(44) 


7(44) 


0-2:7(44) 


0-2,7(44) 


7(44) 


This completes the proof of Theorem 
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5. The heat content of a tt-sector in a ^-wedge. 

In this section, we compute the heat content of a ir-sector in a ^-wedge which share 
one common edge (and vertex). This is a crucial ingredient in the computation of the 
heat content of the fractal polyhedron Ds , which was constructed in Section and will 
be used in Section!^ 

Lemma 11. Let B.,y[R) C such that B-n{R) and share one common edge (and 

vertex). Then, for t ^ 0, 

(45) / dx dyp{x,y,t) 

2 

= - ^t^^^ + {4:'Kt)~^^^ f d-i/i (arcsini/) — i/)) [ drr^e~^ 

2 VTT Jo Jr 

We remark that the coefficient of t is equal to 0 in this case. 


Proof. Similarly to the proof of Lemma the left-hand side of (451 equals 

r-iT ^ 

(47rf)“^ / dOr / dd 2 / dri / dr 2 (rir 2 )e 

Jo Jo Jo Jo 

p'K nl'K pR POO 

= / dOi / d92 / dr\ / dr 2 {rir 2 )e 

Jo Jo Jo Jo 


PIT pItT pR poo 2 

— (47rf)“^ / dOi / d92 / dri / cir2(rir2)e“^''i 

Jo J^ Jo Jo 
(46) ^■.Mi+M2. 

Now 


-(i-I+’'2)/(40 + 2r-ir'2A/(4t) 


PTT PATZ pR poo 

Ml = (47rf)“'^ / dOi / d 62 / dri / dr2{rir2) 

Jo Jo Jo Jo 


,-(rf+r|)/{4t)+2rir2A/(4t) 


(47) = 

and letting r 2 — riA = p, we have that 

p27r pR 


PTT pATT pR poo 

M 2 = — (47rf)“^ / dOi / d92 / dri / dr 2 {rir 2 ) 

Jo J^ Jo Jo 

P7T pO pR poo 

— —(47rt)~^ / d6\ I d 62 I dri I 

^0 J-^ Jo Jo 

PIT pO pR poo 

~—{47vt)~^ / d9i / d92 / rdr dp{p-\-Ar)f. 

Jo J-^ Jo J-Ar 


= _Ji(l_e-«V(40) 


(48) 

We also have that 


PTT pO pR poo 

(47rt)~^ / d9i / d92 / dr j dpAr 

Jo 4-A Jo J-Ar 


-p2/(4t)-r2{l-A2)/(4t) 


2g-pV(4t)-r-‘(l-A^)/(41)^ 


PIT pO pR poo 

{4:7vt)~^ / d9i / d92 / dr dpAr 

Jo J-^ Jo J-Ar 

PIT pO pR poo 

- —{47vt)~^ / d9i / d92 / dr dpAr 

Jo J-^ Jo Jo 

/ TT ^0 pR pr 

d9i / d92 / dr dp A 


2g-pT(«)-r-^(l-A^)/(4t) 


2g-pV{4t)-r^(l-A")/(4t) 


2^2g-A2p2/(4t)_^2(i_^2)/(4j) 
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(4^t)-V2 

pTT 

1 dOi 

r d92 

/' 

drAr^e ^ 

2 J 

0 

J-i 

Jo 


pTT 

{ 

■0 r 

R 

poo 

- (47rt) / 

dSi / 

d92 / 

dr 

1 dpA^r^t 

Jo 

J- 

-f Jo 


Jo 

pTT 

{ 

■0 r 

R 

poo 

+ (47rt) ^ / 

dSi / 

d92 / 

dr 

j dpA^r^t 

Jo 

J- 

-f Jo 


J r 

(4^t)-i/2 

pTT 

1 dOi 

r d92 

f 

dr {A + A|)r 

2 J 

0 




f-TT 

[ 

■0 r 

R 


+ (47rt)“^ / 

dd^ / 

d92 / 

dr 

^ r. 

1 dpA'^r'^t 

Jo 

J- 

-f Jo 


J r 


2 2-A^p^/(4t)-r^(l-A^)/(4t) 


2 2-A^p^/(4t)-r^(l-A^)/(4t) 


^•^^-r\4-A^)l{4t) 


2 2-A^p^/(4t)-r^(l-A-‘)/(4t) 


(49) =: A^i + A^2- 

Now Ni equals 
(47rt)"^/^ 


f dOi [ dOi [ dr (cos(6li - 62 ) + \ cos(ei - 

^0 J-^- Jo 


(47rt)"^/^ 

2 

(47rt)"^/2 


(47rt) 


(47rt) 


(50) 


4 . 

■v/TT 


p7T 

cf rH 

/ d9^ 1 

' d92 / 

dr (cos(6'i + 6^2) + 1 cos(0i 

Jo J 

0 ^0 


pTT 

/■f+Si 


/ d9^ 1 

1 dp 

/ dr (cosp + \ cos p\)r‘‘e 

Jo J 

e^ 'J 

’0 

r d9^ 1 

f" dp r 

J 2 —(sin ?7)^/(4t) 

dr r cos 77 e ^ u / y j 

Jo J 

61 Jo 



/■f+®1 

fH 2 

- / d9i 

/ dp 

/ dr (cost; + 1 cos77|)r e 

Jo 


Jo 

1 pTT 

cf+01 

fR 2 

T; “'J 

1 dp 

9i 

/ dr (cos 77 + 1 cos77|)r e 
/o 

f" d9 f 

— pR 

^ T / 1 2 —(sin r7)^/(4t) 

dp drr cos77e ' u ; 

Jo Je 

Jo 



pi 

poo 

' y. 

^^ 7 /J (arcsinT/j — t/j) / c/rr' 

Jr 


N 2g-T-'^{sin(6li+e2)) /{4t) 


2 -r (aiup) /(4t) 


(ainij) /{4t) 


+ UO, 

by integrating by parts with respect to 6 . In addition, as for the computation of I 5 (see 


(33 1 ), we have that 

r7T rO pR poo 

N 2 = (47rt)~^ / d 6 i / d 62 / dr dp A 

Jo J-^ Jo Jr 


2^2g-A2p2/(4t)-r2(l-A2)/(4t) 


2 t 

TT 

t 

TV 

t 

TV 

t 

TV 

t 

TV 


pir pO poo 

J d9i J d 62 j dp 

p7V pO 

I */-, 

f de^ r 

Jo Jo 


+ 0(16-”"/^"^*’) 


(A2p2 + 1-A2)2 

d 62 ( —(cot(0i — 62 ))^ + arctan(tan(6>i — @ 2 ))^ + 0 {te~^ 

(sin(^ 6 /i — 6 ^ 2 )) / 

2 , cos(Si + ^ 2 ) a,rctan(tan(0i + 02 ))^ + 0(te~^ 


d02(-(cot(0i+e2)) 


[ d9i [ djy r—(cotr;)^ + arctan(tan» 7 )^ + 0(te ^ 

Jo Je, \ (sinry)3 J 

f dO f dp (—(cot arctsin(ta.np)\ + 0 {te~^ t ^ 0 . 

Jo Je V (sm?7)4 J 
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Note that 


arctan(tanJ 7 ) = -f 

ri-n if?7G(2,V)- 


Hence 


N 2 = — f dd f dr] [—{cotrj)^ + arctan(tan 77 )^ + 0{te ^ 

Jo Jb \ (sin 77 )'^ J 


f^de f"drji-{cotvr + ^"''""^ 

Jo Je 



t 

f 2 

/■f+« / 

+ 


/ dO 



TV , 

JO 

4 V 


t 

pTT 

/■f+® / 

+ 

/ de 

/ f ■ 


TV , 

h ■> 

'e V 


(siniy)^ 

2 (>7 - tt) cos rj 

(sin 77)3 

si2 (77-71) cos 77 


(sin 77)3 

+ ^ + - 

2tt tt \ 4 J 


+ 0{te 


-RVm-s 


(51) 


t t , t / n 

71 


= ^t + 0{te ^ tiO. 


Thus, by (481, (491, (501 and (511, we have that, as t ), 0, 
(52) 

M2 = -^^' 

VTT 


+ (47rt) f (i'0 (arcsini/i —'0) f drr^e ^ ^ + 0(te 

Jo Jr 


■r 2 cos 92r/(4:t) — {ri sin0i—r 2 sin 02)^/(4t) 


Therefore, by (461, (471 and (52l, astJ^O, 

p-n /’^ f°° 

(47it)“^ / ci0i / d6l2 / dri / cir2(rir2)e“^’’^ 

Jo Jo Jo Jo 
= -+ (47it)~^^^ f d?/) (arcsin 7 /) — 7 /)) f drr^e~^ + 0{te~^ /(^*)y 

~ Jo Jr 




□ 


We note that 


V^{t-,R)=j dx dyp{x,y,t) 

^ Jb^(r) Jw^ 

2 2 

= / dx dyp{x,y,t)+ / dx dyp{x,y,t). 

J B^{R) JBn(R) 

2 22 

By Theorem and Lemma [m this implies that, as t J 0, 

/ dx dyp{x,y,t) 

Jb^{r) Jw^ 

= ^-Tii/2 + 

4 0F V 2 J 

+ (47it)~^^^ f dtp {arcsinip — Ip) f drr^e~^ + 0{te~^ /(8‘)^ 

Jo Ji? 


(53) 


4 x/Tr TT 


+ (47it) f dtp {arcsintp — tp) f drr^e ’’ + 0(te 

Jo Jr 

The result of Lemma 11 and formula (531 will be used in Section]^ below. 
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6 . The heat content of the fractal polyhedron Ds. 

In this section, we use Theorem[^to compute the heat content of the fractal polyhedron 
Ds which was constructed in Section To do this, we adapt the scheme of [ 6 ] to the 
three-dimensional setting below. The key step in [ 6 ] was to obtain a renewal equation by 
making a suitable Ansatz for the heat content. The corresponding Ansatz has been made 
here in (601 and (611 for 0 < s < \/2 — 1 , s 7 ^ |, and in (631 and (651 for s = | • 

In order to derive the required renewal equation, we need to compute the contribution to 
the heat content Hds (t) from Qo and Qi.i- We do this below in Lemma [T^ and Lem ma|13| 
respectively. In what follows, for A C d{x,A) is the 3-dimensional analogue of (181. 
We make the following approximations for UDs{x;t). 

Let 0 < (5 < min | ^| and let x € Ds- If d{x, dDs) > 3, then we have that 

|MD,(a:;t) - 1| < 

by the principle of not feeling the boundary, O Proposition 9(i)]. We define 
F = {x G Ds : d{x,dDs) < 5,d{x,e) > 5 for all edges e G dDs}. 

If X G F, then we have that 
( 54 ) \uDsix-,t) — UH{x-,t)\ < 

where 


Unix'jt) = (dirt) / d<}e 

J — d(x,dDs ) 

i.e. 7L is a half-space whose boundary contains the face of dDs nearest to x. Let 
E = {x G Ds : d{x, e) < 3 for some edge e G dDs, d{x, v) > 3 for all vertices v G dDs}. 

If X G E, then we have that 

(55) \uDs{x\t) — uw{x',t)\ /(s*)^ 

where W is the infinite wedge W:k for entrant edges and W is the infinite wedge Wsji for 
re-entrant edges. (See the proof of Lemma 12 for further details). The estimates (541, 
(55 I follow by similar arguments to those given in the proof of Lemma with D = Ds, 
E = H, W, E = F, E respectively and G = {x G Ds : d{x, E) < d}. 

It remains to approximate UOs{x;t) for x near a vertex of dQo H dDs, 9(5i,i H dDs 
respectively. We only require the contribution to the heat content Hosit) from these 
vertices to derive the required renewal equation. The relevant approximation to make 
here is via a one-sided infinite cone with vertex v G dDs such that dCv ^ {x G dDs : 
d{x,v) < (5}. Definition aside, no viable expressions are known for ucs{x;t) in this 3- 
dimensional setting. For our purposes, it is sufficient to approximate the neighbourhood 
of each vertex v G {dQo U clQi.i) n dDs by a cube S'„. Each cube S„ centred at v has 
side-length 23 and is chosen such that the faces of dSv are pairwise-parallel to those of 
dQo. We are interested in the contribution to the heat content Flosit) from the region 
Sv C] Ds. There are two cases to consider. Either the vertex v is entrant and 5'„ n Ds is 
I of Sv, or the vertex v is re-entrant and Sv n Ds is | of 5„. If u is entrant, then the 
coefficient of in the expansion for Hds (t) is equal to -Gr /2 by separation of variables. 
Unfortunately, we were unable to compute the coefficient of for a re-entrant vertex. 
However, the contribution to the heat content (t) from each region Sv n Ds is of order: 
3^ from the volume, from the surface area of the adjoining faces of the vertex, and 

3t from the length of the adjoining edges of the vertex. Thus, if we choose 3 as follows: 

(56) <5 = 8t'/2(log(t-i))i/", 

then the contribution to the heat content Hosit) from each region Sv H Ds is 
0{d^^{Iog{t~^))^^^). We note that this choice of 5 gives = 0(t®). 

Lemma 12. Let 0 < s < y/2 — 1. Then 


1/2 
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Proof. Partition Qo into the following sets. 

(i) dQo n dDa has 32 vertices; Vi,i = 1, • • •, 32. At each vertex Vi, consider a cube Si of 
side-length 25 centred at Vi. Let St = SiD Qo, i = 1, ■ ■ ■ ,32 and S = 

(ii) dQo n dDs has 36 edges; ej,j = 1,..., 36. Let Ej = {x G Qo ■ d{x, Cj) < S,x ^ S}. 

(hi) F= I® € Qo : d{x,dQondDs) < 5,x ^ 

(iv) The interior of Qo minus (i), (ii) and (hi); an open polygon P^ with distance at least 
5 to dDa- 

(v) The remainder, which has measure zero. 

The contribution to the heat content from (iv) is \Ps\+0{e~^ /(8‘)^ = |p^|0. 

To compute the contribution from (ii), there are two types of edges to consider. Each 
edge Cj is the intersection of two faces of dDs. For fixed j, let IIj denote the plane 
which is orthogonal to these faces and intersects Cj in exactly one point. Apply Lemma[^ 
with D = Ds n IIj , F = respectively, E = Ej C] flj and G = {x G Ds n Ilj : 

d{x, Ej n flj ) < 5} so that either; 

(I) the contribution from Ej n IIj can be approximated by that from a sector Bs{^) in a 
wedge VF:|, or 

(II) the contribution from Ej n Ilj can be approximated by that from a sector Bs (tt) in a 

wedge . 

2 

We can use Lemma 1^ to deduce that the contribution from the edges of type (I) is 

12(1 - 25) lEj n n, | - 25— + -t] + 

\ VtT TT y 

= 12 |.B,| - 25(1 - 2(5)^-h-(1 - 2<5)t + t fH. 

\ VtT TT y 


In addition. Lemma 11 gives that the contribution from the edges of type (II) is 

24(s - 25) i^\E, n n,| - 5 —j + 0(t®/^) 

= 24 iPEj] - 5{s - 2-5)^ j + 0(t®/^), t i 0. 

Each cross-section of E is a union of rectangles and cuspidal regions. Thus by Section]^ 
we have that the contribution from F equals 


|E| -U\dF) 




uo. 


For edges of type (I), the sector Ej n Ilj has two cuspidal neighbours. For edges of type 
(II), the sector Ej nllj has one cuspidal neighbour but by Lemmajllj the term of order 
is half that of LemmaEven though the sector and cusp terms of order 1?^'^ cancel out, 
the remainder is dominated by the contribution from (i), which is 0(t®^^(log(t“^))®/^). 
This completes the proof of Lemma □ 


Following the strategy of [B], in order to compute dx uosi^',!), we introduce a 

model solution which approximates ud, in one of the six components of Ds — Qo. Consider 
the half-space H = {{xi,X 2 ,X 3 ) G R® : < 0} and attach one of the six components of 

Ds — Qo to H. The resulting set is 


Hs = interior < 

HU 

U U 


[ 

J>ll<i<^N(j) 


Let uhs be the unique solution of ||^ and ([^ with D = B/s. Define 



dxuHs 
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Applying Lemmaj^with D — Da, F — Ha, E — Ha — H and G = {x G Ha : d{x, E) < e}, 
where e = (1 — 2s — s^), we have 


(57) 



—Qo 


dx ud 



dx un, [x] t) 
dxuHa{x\t) + 0{e~‘‘ 


= 6£(t) + UO. 


In contrast to [^, where the temperature of the boundary is fixed for all f > 0, we must 
account for the fact that Ha — H and Ha — Da feel each other’s presence. The choice of t 
above is a lower bound for the distance between Ha — H and Ha — Da- 
Similarly to Lemma we have; 

Lemma 13. Let 0 < s < y/2 — 1. Then 

f 12(5 

/ dxuHa{x-,t) = — 5s^(l — s^)—= H- 1 + h{t), t 4. 0, 

jQi.i 

where \h{t)\ < for t < and some constant C > 0. 


Proof. The proof of Lemma[^is analogous to that of LemmaWe note that in this case 
there is an additional type of edge to consider. Namely, the edges where the contribution 
from {E n Qi,i) n Dj can be approximated by that from a sector (5) in a wedge . 
There are 4 such edges so the contribution to the heat content is 


+1/2 1 \ 

4(s - 25) ( KS n Qi.i) n HjI - 5— + -t)+ 

' y/TT TV J 

= 4 (|£; n Qi.il - &{s - 25)^ + + 0(t3/"), t 0, 


by (53l. This completes the proof of Lemma 13 by our choice of 5, (|56|). 


□ 


Below we state and prove the corresponding 3-dimensional result to [6l Proposition 4] 
for completeness. 

Lemma 14. Fix 0 < s < y/2 — 1. Then 

(58) E{t) = 5s^E ( -I- — 5s^(l — s^)^F^ _|_ -|- h{t), 

J TT 

where \h{t)\ < C't®^^(log(t“^))®^^ for t < s^ and some constant C > 0. 


Proof. Re-write E{t) as 
E{t) = 


I 

Qi,i 


dXUHa {x',t) + 


L 


Ha-{HUQi,x) 


dx uh, (x; t). 


Then Ha — {H U Qi.i) consists of 5 copies of Ha — H scaled by a factor s, say Ai,..., A 5 . 
Each of these copies has a face fi connecting it to Qi,i. Let Ha^ be the half-space such 
that OHai 3 fi and J/a; 3 Qi,i. Put Fi = AiU HAi U fi. Then Fi is a copy of Ha scaled 
by a factor s. Thus, by scaling, we have 

/ dxuFi{x-,t) = s^ / dxuHa{x\t/s^) = s^Eft/s^). 

J Ai J Ha-H 

Define Gi = {x G Ha : d{x,Ai) < e}, where e = 27 i-a)’ ^ Lemma with 

D = Ha, F = Fi, E — Ai and G = Gi, we have 

/ dxuHa{x;t)= / dxuFi{x-,t) + 0{e~'^ t fO. 

J Ai J Ai 
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Hence 


L 




dxUH,{x-,t) = ^ / dxUH,{x\t) 


L 

5 f dxuFi{x-,t) + 0{e~'^ /(8*)^ 
J Ai 


= 5s^Bit/s'^) + t i 0. 

Combining this with Lemma[^ gives the result. 

We must now consider the different regimes for s. 


□ 


Lemma 15. Let d = | + 


_ 3 j_ log 5 


2^2 log s 


and fix 0 < s < %/2 — l,s 7 ^ 4. Then there exists a 


periodic, continuous function ps : R —>■ R with period log(s such that as t fO, 
s^ 5s^(l - s^) 12s 


(59) E{t) = 


1 — 5s^ 1 — 5s2 ^/tv 7r(l — 5s) 


t + ps{logt)t‘^ + 0(t®/^(log(t ^))®^^). 


Proof. Define 
(60) 


qs{t)t‘^ = E{t) - 


s“ . 5s^(l - s^) 12s 


+ 


1 — 5s® 1 — 5s2 '^(1 “ 5s) 

Substitute (581 into ( |60[ ) to obtain 

1 - s^ \ fi/® 


qa{t)t = 5s E ( ^ ) — 5s 


1 — 5s® 


+ 25s* 


— 5s 


1 —5s®y -^tF \7r(l — 5 s) 


12s 


t + h{t). 


From (601, we also have 
t \ t‘^ 


s® / 5s® 


= E 


s^ 5s(l-s^)t^/^ 12 

® ^ 1 — 5s® 1 — 5s® y^r 7rs(l — 5s) ’ 


which implies 

or equivalently, 


H = <lsit)t - h{t), 


qs{t) = Is ( ^ ) + ■ 


Similarly to [51 Proposition 5], define Ps (log t) by 


(61) 

Then 

hence 


qs(t) = Ps(logt) - ^/i(ts^®)(fs^®) 


j=i 


Ps flog 4 ) = Psit) + '^h(ts^^){ts'^^y‘^ = ps{logt). 

V ® / 7=1 

Since there is a constant C > 0 such that \h{t)\ < C't®/®(log(t~^))®/^ for t < s®, there is a 
constant (7 > 0 such that 

00 / 00 00 

^/i(ts"®)(ts"®)-‘' < C't(®/®)-‘* (log(t-'))®/"^(s®-''‘*)^ + (-logs)^/"^(2jf/"(s®-"'")® 


j=i 

(62) 


1=1 


1=1 


= 0(f{®/2)-d(log(t-l))®/2). 


Combining (601 with (611 and (621 gives (|59|. 

Lemma 16. Let s = |. Then there exists a periodic, continuous function pi : R - 
with period log 25 such that as f ), 0, 

1 A 19 


□ 
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Proof. For s = |, d = 1. Define 
(63) 


, , , 1 6 6 , 12 

qi (t)t = Bit) -^ - — H -^ 1 log t - 1 . 

5^ 120 25 57rlog5 Stt 


Substitute (581 with s = | into (631 to obtain 

1 16 6 
q.it)t = -Ei25t) + g^tlogt + h(t). 


By considering (631 with t replaced by 25t, we obtain that 

qi (25t)t = qi (t)t — h(t), 

5 5 

or equivalently, 

(64) 

Dehne pi (logt) by 


qi{t) = qi{25t) + h{t)t 


(65) 


qi (t) =pi(logt) - ^h(t(25) ■’)(t(25) \ 

5 5 ' ^ 

j=i 

we have pi (logt) = pi (log25t) via Since 


15 


Then, as in the proof of Lemma 

holds for s = |, we obtain the same remainder estimate as in Lemma 


15 


□ 

Combining Lemma |12[ ( |57[ ) and Lemma ] 15[ Lemma|16|we obtain ( |16| l, ( |17[ ) respectively. 
Using the fact that t i—> Hosit) is continuous, it can be shown that t i—>■ p 3 (logt) is 
continuous for 0 < s < y/2 — 1, see [6l Section 5.1]. 
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